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Gauge field theory of transport and magnetic relaxation in underdoped cuprates 
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Based on recently proposed (7(1) x 517(2) Chern-Simons gauge field theory, an interpretation of 
the transport and magnetic relaxation properties of underdoped cuprates is proposed, taking into 
account the short range antiferromagnetic order. The interplay of the doping-dependent spin gap 
(explicitly derived by us) effect and dissipation due to gauge fluctuations gives rise to a crossover 
from metallic to insulating behavior of conductivity as temperature decreases, in semi-quantitative 
agreement with experimental data. For the same reason the magnetic relaxation rate shows a 
maximum nearby. Various crossover temperatures related to spin gap effects are shown to be 
different manifestations of the same energy scale. 



PACS Numbers: 71.10.Pm, 11.15.-q, 71.27.+a 

The understanding of the anomalous normal state 
properties of oxide superconductors has been a chal- 
lenge for theorists since their discovery. Recently 
great attention ias been focused on underdoped 
superconductors, uu where the pseudogap (spin gap) 
effects are essential. We will concentrate on the dop- 
ing range where the short range antiferromagnetic order 
(SRAFO) exists, B and propose an interpretation of the 
transport and magnetic relaxation properties in this re- 
gion, based on the recently proposed C7(l) x SU{2) gauge 
field theory O 

The linear temperature dependence of resistivity in 
most of oxide superconductors over a wide range of tem- 
perature is well establishecLand a number of explana- 
tions have been proposed, B including the U{\) gauge 
field theory. H On the other hand, in underdoped sam- 
ples, a resistivity minimum and a crossover froip_,rpntal- 
lic to insulating behavior has been observed. E3"E3 A 
similar divergence of resistivity at low temperatures has 
been found ixL,superconducting samples in strong mag- 
netic fields, 113 suppressing superconductivity. An ap- 
parently "obvious" explanation of these two related phe- 
nomena would be localization of charge carriers in two 
dimensions. However, a morcLiiarcful comparison of the- 
ory with experiments shows O that the localization ef- 
fects including carrier interactions cannot correctly inter- 
pret the data. Several other explanations have been pro- 
posed based |On non- Fermi liquid (FL) behaviowif-charge 
carriers, [EjO but the zero field experiments llJ^llJ have 
not been addressed, except for [|^ where a gauge field 
approach has been used. We, instead, will concentrate 
on the latter case. We will show that the presence of 
SRAFO, leading to a finite mass of spinous (bosons) 
is the correct starting point in this doping range. The 
self-generated U(l) holon-spinon [h/s) gauge field bj 
comes singular due to coupling with holons (fermions) 
which, in turn, renormalizes the massive spinous in a 
nontrivial way. At low temperatures, effects due to fi- 
nite spinon mass prevail leading to insulating behavior. 



while at higher temperatures the dissipation caused by 
the gauge field dominates and gives rise to metallic be- 
havior. For similar reasons, the spin relaxation rate is 
low at both low and high temperatures, reaching a max- 
imum near the resistivity crossover point which is also 
consistent with experiment. H 

Following a strategy previously applied to the ID i— J 
model which hasjusproduced there the known exact Bethe 
Ansatz results, O the Cjaern-Simons bosonization with 
U{1) X SU{2) gauge field H was applied to the 2-1) t- J 
model in the limit t ^ J, allowing us to rewrite the par- 
tition function (and the correlation functions) in terms 
of a spin i fermion field ipa, ct — 1,2, minimally coupled 
to a U{1) field B (gauging global charge), and an SU{2) 
field V (gauging global spin) whose dynamics is given 
by a Chern-Simons (C.S.) action. We decomposed the 
fermion field '0a h^^to product of a spinless fermion field H 
(holon) and a spin 2 boson field S^ (spinon), satisfying 
the constraint S^Sq = 1, thus introducing a local f/(l) 
gauge invariance called h/s. We proved the existence of 
an upper bound of the partition function for holons in 
a spinon background, and we found the optimal spinon 
configuration (s-|-id-like RVB state) saturating the upper 
bound on average. After neglecting the feedback of holon 
fiuctuations on field B and spinon fluctuations on field 
V, the holon field is a fermion and the spinon field is a 
hard-core boson. Within this approximation the "mean 
field" (MF) B produces a tt flux phase for holons, con- 
verting them into Dirac-like fermions, while the V fleld, 
taking into account the feedback of holons produces a 
gap for spinous vanishing in the zero doping limit. 

The continuum action for AF fluctuations around the 
"MF", described by a spin ^ boson field Za,a = 1,2 ( 
still "spinous" ) is given by: 



dxod^x[v;^\ido - Ao)z\^ -\{d^- A^)zf^ 



(1) 



where A is the h/s gauge field, g = 8/ J, Vg — ^/2Ja, 
with a the lattice constant. The spinon "mass" term 
rn^ ~ (y^) ^ —d\nd (the main novelty) is due to aver- 
aged perturbation caused by holons of concentration S via 
V. This explicit doping dependence was derived, rather 
than assumed in the theory. It produces a SRAFO, with 
correlation length ^af ~ (^'^1'^4)t"^' f^l^Y consistent 
with the neutron scattering data. O 

Neglecting the gauge fluctuations, holons are described 
by FL theory with a Fermi surface (FS) consisting of 4 
"half-pockets" centered at (±|,±|). The MF B turns 
the spinless fermion H into two species of 2-component 
Dirac fermions ip^^'jr = 1,2, each of them being sup- 
ported on one Neel sublattice. The continuum action for 
these fermions is given by: 

/dxod'x^ V5('-)[7°(9o - erAo - S) + t{^ - er4)]^^"\ 

-' r 

(2) 

where 4 = 7m^M'^ = 7m^M'7o = 0-2,7^ = {ay,aa:), the 
charges e^ = ±1, depending on sublattice. After inte- 
grating out the gapful Dirac modes, we end up with FL- 
like system of holons with Fermi energy ep ^ tS, inter- 
acting through gauge field A. 

As shown in [Q , for the gauge field model the in-plane 
resistivity is approximately given by 

R = limc.[(Imn^(c.))-i + (Imn^(c.))-i], (3) 

where 11^ and 11^ denote the transverse polarization 
bubbles (at g = 0) due to the h/s currents of holons 
and spinous, renormalized by gauge fluctuations. The A 
propagatoE-for small \q\, uj, uj/\q\ in the Coulomb gauge is 
given by: iO 



{A^A^){q,L0)^{tiuX,iq)+xm 
{AoAo){q,uj) ^ {vh +cjp)~\ 



2\-l 



(4) 



where A^ ^ K,/\q\,K ^ 0((5) is the Landau damping due 
to a finite FS for holons, x ^ Xh + Xs,Xh ^ mj^^ ^ 
0{5~^),Xs — VsTTi^^ ^ 0{{—dlu6)~^) is the diamag- 
netic susceptibility, i^h is the holon density at the FS and 
ujp is the plasmon gap. 

An estimate of the holon contribution to resistivity can 
be derived as in B, 



Rh ~ {eFn„rpr' + {-f'V^F, 

X 



(5) 



where Tj,„p is the transport relaxation time due to impu- 
rities. 

To estimate the spinon contribution, we derive the 
large scale behavior of the spinon current 7'^' = z* D^z 
correlation functiDu, where D^ = dfj, — A^, by eikonal 
approximation, ^3 strictly preserving gauge invariance. 



We use spinon Green functions at zero temperature, as 
partially justified by the spinon gap, but we retain the 
temperature dependence of gauge fiuctuations. 

We apply the Fradkin representation cil to the spinon 
propagator {z{x)z*{y)) — G{x,y\A). It can be de- 
rived using a first-quantized path integral form of the 

propagator, with metric ( h +), replacing integration 

over trajectories q^{t) by integration over 3-velocities 
4>ti, — 9^(0' M = 0: I7 2. Rescaling .tq to VgXQ one obtains: 



G{x,y\A)^i 



-ism \„is{dfj,-Af^) 



{x,y) 



I?(/.^(t)e*/o°*'(*)'^* 



^ ^ij/;^A^(t)<t,'^(t)dt I" ^3 ip^ix'^-y'--J^4,t-(t)dt) 



(6) 



Using an identity (Eq. (41) in the second paper of p4|), 
the integral /^ A^(t>^'{t)dt with i^ = A^{x + J* (j){t')dt') 
can be decomposed into a sum of an integral along a 
straight line (denoted by J^ A) and a gauge invariant part 
depending on the field strength F^i,. Thus G{x,y\A) = 
exp{— i/^ A}G(x,y|F). The spinon current density cor- 
relation {j^(x)j'^{y)), /i = 1,2, is approximately given 
by {D'X{x)G{x,y\A)D'X{y)G{x,y\ - A)), where (•) de- 
notes average w.r.t. A. The gauge-dependent terms of 
the two spinon propagators exactly cancel each other, 
yielding a strictly gauge-invariant result, at large scale 
given approximately by {■£-G{x,y\F)-^G{x,y\ - F)). 
The A-^- average involves contributions weighted by 
"magnetic field" correlations {Ff^u{z)Fpa(w)) , ji, v, p,a — 
1,2, approximately evaluated for jz" — w°| ^ T~^ 
as in gH, obtaining {5^p5,, - J^,^^^^)^e-I?--I^.^g2^ 

where q^ = {-^)^ is a momentum cutoff related to 
the anomalous skin effect due to the Reizer singular- 
ity in the A-^ propagator, o The Aq average involves 
contributions weighted by "electric-field" correlations 
{Fq^[z)F^i,{'w)Y fi,v = 1, 2. Since they vanish in the limit 
q,uj ^ (see (Q)), their contributions will be neglected. 

Integrals in (|g) can be approximately calculated for 
relatively low temperatures (T < X"^s) ^^'^ ^^^ current- 
current correlation becomes: 

CfX^ 

_r?|g(i72|S|TO2/^2_|-|2> n2\-l/2n2 /^s 

.e ix ^ ^=^0 1^1 ^(a;^ - Ixp) ^'^\ (7) 

where, for a real argument, / is monotonically increas- 
ing, vanishing at zero argument and g is monotomically 
decreasing, vanishing at large arguments. Their expHil 
expressions are lengthy and will be given elsewhere. O 
In deriving the spinon current correlation at g = we 



evaluate the x-integration by saddle point. For Xq ^ q^^ 
the integral is dominated by a complex saddle point at 
\x\ = 2q^ o:{xo), with finite a{xo) (in the first quadrant), 
having a weak dependence on xq- To justify the saddle 
point approximation we need to assume T > xiTT-sQa- It 
turns out that in the physical range of parameters con- 
sidered in the paper, this and the above conditions are 
both satisfied for temperatures between tens and a few 
hundred degrees. 

Let us define 11+ (w) ~ J^ 



'dxoij^'jTiq 



0,xo) 



\ Using the Lehmann representation we find 



lim^^o Imn-^(w)w"i = -2£jRen+(0). Taking note that 
the main contribution comes from small a;*^, introducing 
lower cut-off and performing scale renormalization, we 
obtain for the u ^ Q limit: dKelV^ {lo) / doj 



Re 



X 



where 



\Z\e- 



^fiam, r^fiam. (s) 

A 



(renormalization eliminates the contribution of the g 
function, being subleading). We find that at large xq 
arga(O) — j and arg/"(Q;(0)) = 0. We extrapolate 
a{xo) to a(0), keeping these features. This way we re- 
cover the correct behavior, R -^ oo, as T — > 0. As 
xo — > 0, the saddle point extrapolates to Xs ^ q^ e*~, 
and we find the "spinon contribution" to resistivity: 






|a(0)r 



, \Z\"' 
sin(6'/4) ■ 



(9) 



In Fig. n our calculated resisitivity (sum of (pi) and 
(y)), is plotted as a function of temperature for various 
dopingsjji comparison with experimental data taken on 
LSCO M (inset). We have taken t/J = 3, J = 0.1 eV. 
Apart from the resistivity scale, there are no other ad- 
justable parameters (similarly for Fig. g). We find a 
resistivity minimum below lOOK in very good agreement 
with experiment. We see from (H) that the imaginary 
part of Z is proportional to temperature T . At low tem- 
peratures the spin gap effect {'^ rus) dominates, — > 0, 
so the system shows an insulating behavior. (The func- 
tional dependence R ^ ^/T, different from the "stan- 
dard" exponential law due to spin gap, is a prediction 
of our theory.) On the contrary, at higher temperatures 
the imaginary and real parts of Z become comparable, 
so the resistivity grows with temperature due to gauge 
fiuctuations via \Z\. Moreover, the minimum shifts to 
higher temperatures, as the doping decreases, also in 
agreement with experiment ( our theoretical prediction 
ml ~ —(Sin (5, rather than ~ (S is responsible for this 
shift). We have also compared the calculated conductiv- 
ity in the semi-log scale with data taken oiua^very good 
single crystal of Lai.96Sro.o4Cu04 (inset). B'EJ We find 



a symmetric shape of curve around the maximum, and 
a reflection point as well as a linear piece on the low 
temperature side in both theory and experiment. So far 
we have not included the external magnetic field. We be- 
lieve the experimentally observed crossover from metallic 
to insulating behavior in strong_magnetic fields when su- 
perconductivity is suppressed, 113 can be understood in a 
similar way, and this issue will be addressed in our future 
communication. O 
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FIG. 1. The calculated temperature dependence of in-plane 
resistivity (sum of (H) and (H)) for various dopings S in com- 
parision with the corresponding experimental data (inset) on 
La2_5SriCu04 in units of mQcm, taken from ]10| . 

Now turn to the spin-lattice relaxatipn rate Tj^ which 
can be expressed approximately as: 119 



(TiT)- 



lim / (fqT{q} 



lmxsiq,uj) 



where Xs is the spin susceptibility and J-{q) is the form 
factor. To evaluate Xs we use the representation for spin 
deduced at large scales S^ ~ e^'^^^^z*az{x){l — ph{x)), 
where p^ is the holon density, to be replaced by its aver- 
age S. Around the AF wave vector Qaf = (tt, tt) we find 
that {S{x) ■ S{0)) ~ (1 - <5)2e"l^l (G(x, 0|F)G(x, 0| - F)) 
which can be calculated as before. Define x+(q', w) = 
/p dxo{S ■ S)*{q,xo)e^'^"'^, using the Lehman represen- 
tation and taking into account that J^{q) is even in q 

one obtains (TiT)-^ = -2/d^g.F(g) ^^°^t,<5-"=°^ Since 

J-{q) is peaked around Qaf for Cu, integrating over q 
in a small region around that point, and introducing a 
cutoff in the real space A '^ t:/\xs\, we find 

{T,T)-' ^ (I - S)'V~S\Zri{acos{^-) + 6sin(^)), (10) 

a = ReX,dVo(2|y|a(0)),& = -lmJ^d^yM2\y\aiO)), 
and Jo is the zero-order Bessel function. 
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imum of the in-plane resistivity), the spin gap crossover 
temperatures, detected by NMR Tq and by resistivity 
T*, identified with their respective reflection points, in 
the low-doping region ( S ~ 0.02 — 0.08). The last two 
temperatures, roughly speaking, limit from the above the 
region of significant spin gap effects and validity of our 
approximation. These crossover temperatures are differ- 
ent manifestations of the same energy scale. The fact 
that their relative order T* > Tq > Tm-i, as well as 
their order of magnitude agrees with experiments, pro- 
vides further support for our approach. 

The work of P.M. was partially supported by TMR 
Programme ERBFMRX-CT96-0045. 
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FIG. 2. The temperature dependence of calculated 
spin-lattice relaxation rate {TiT)~^ given by (|lCl|). Inset: 
•^^(TiT)"^ in the Cu02 planes of YBa2Cu306.52 single crys- 
tals in units of s~^K~^ , taken from llSl. 

In Fig. H we plot our calculated spin-lattice relaxation 
rate (TiT)^^ for ^^Cu as a function of temperature for 
various dopings in comparison with expetimental data 
taken on underdoped samples of YBCO. M We observe 
a maximum near the crossover temperature for conduc- 
tivity, although the shape around maximum is not sym- 
metric anymore, due to the presence of the cosine term. 
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FIG. 3. The calculated metal-insulator crossover tempera- 
ture Tm-i, the reflection point Tq in NMR (rir)"\ and the 
reflection point T* of R as functions of doping. 

To summarize, we have shown using the C/(l) x SU{2) 
gauge field theory that the metal-insulator behavior 
crossover and peculiar behavior of NMR relaxation in un- 
derdoped cuprates might be due to the interplay of the 
spin gap (derived in our approach) effect and the gauge 
field fluctuations. More precisely, the crossover is taking 
place when the real and imaginary parts of Z (Eq. dq)) 
become comparable. In Fig. g we have plotted three 
different crossover temperatures related to the spin gap 
effects, namely the metal-insulator crossover Tm-i (min- 
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